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Abstract 

Coarse-grained theories of dense polymer liquids such as block copolymer melts predict a uni- 
versal dependence of equilibrium properties on a few dimensionless parameters. For symmetric 
diblock copolymer melts, such theories predict a universal dependence on only yN and N, where 
X is an effective interaction parameter, N is a degree of polymerization, and TV is a measure of 
overlap. We test whether simulation results for the structure factor S(q) obtained from several 
different simulation models are consistent with this two-parameter scaling hypothesis. We com- 
pare results from three models: (1) a lattice Monte Carlo model, the bond-fluctuation model, (2) 
a bead-spring model with harsh repulsive interactions, similar to that of Kremer and Grest, and 
(3) a bead-spring model with very soft repulsion between beads, and strongly overlapping beads. 
We compare results from pairs of simulations of different models that have been designed to have 
matched values of N, over a range of values of xN and N, and devise methods to test the scaling 
hypothesis without relying on any prediction for how the phenomenological interaction parameter 
X depends on more microscopic parameters. The results strongly support the scaling hypothesis, 
even for rather short chains, confirming that it is indeed possible to give an accurate universal 
description of simulation models that differ in many details. 



* Corresponding author: 



morse012@umn.edu 



I. INTRODUCTION 



Modern understanding of polymer liquids relies heavily on the study of highly simplified 
models, and of theories that predict highly universal behavior. The best developed example 
of this is the theory of dilute and semidilute polymer solutions in good solvent p~H3]. Exper- 
iments and simulations in this regime are well described by a scaling theory that predicts 
a universal dependence on only two dimensionless parameters, which measure the extent 
of overlap and the strength of the bare repulsion per chain. The theory is based on an 
analysis of a highly simplified model, due to Edwards jl] , of polymers as continuous random 
walks with a short-ranged repulsion between segments. Perhaps the most important step in 
testing this theory was the demonstration that experimental data for polymer solutions is 
consistent with the scaling hypothesis, which was achieved by comparing data from a variety 
of chemical systems, over a range of conditions and chain lengths [6] . 

We are interested here in dense polymer liquids, such as block copolymer melts and 
polymer blends, rather than solutions. Coarse-grained theories of these systems are based 
on a simple generalization of the Edwards model. To describe dense liquids, the model is 
studied in a nearly-incompressible limit in which density fluctuations are strongly suppressed. 
Following Matsen, we will refer to this as the "standard model" [7j for this class of liquids. 
The purpose of the current paper is to verify a scaling hypothesis suggested by a long series 
[51-TT7] of theoretical studies on this model, which predict behavior that depends only on the 
parameters required as inputs to self-consistent field theory and one additional parameter, 
N, which is a measure of overlap. Specifically, we test a prediction that the structure factor 
S(q) in disordered symmetric diblock copolymer melts should exhibit a universal dependence 
on two parameters xN an d N, independent of many details of different simulation models 
and experimental systems. This is tested by comparing simulation results for S(q) obtained 
from several different simulation models over a range of chain lengths and parameter values. 



II. THEORETICAL BACKGROUND 

The vast majority of theoretical work on block copolymers and polymer mixtures has 
thus far been based on the self-consistent field (SCF) approximation [T8l - l2~T] . Self-consistent 
field theory (SCFT) is a coarse-grained density functional theory that expresses the free 
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energy functional of an inhomogeneous polymer liquid as a sum of the free energy of an 
inhomogeneous ideal gas of random-walk polymers and an excess free energy, and in which 
the excess free energy is taken to be a local functional of composition. In the simplest and 
most common form of the theory, for systems with two types of monomers A and B, the 
excess free energy per monomer / ex is taken to be a function / ex = kTx{T)<j)A4>B, in which 
<f>i represents the local average volume fraction of i monomers, and x is t ne Flory-Huggins 
interaction parameter. 

SCF theory predicts a universal dependence of all equilibrium properties on a small set of 
dimensionless parameters. In systems with two types of monomer, such as diblock copolymer 
melts and binary polymer blends, SCF predictions all depend on a product %iV that quanti- 
fies the excess free energy of mixing per chain, where N is a degree of polymerization. Other 
relevant parameters include ratios of chain and block lengths, macroscopic composition (in 
polymer mixtures), and ratios of statistical segment lengths or radii of gyrations. Thus, for 
example, for idealized diblock copolymer melts with monomers of equal statistical segment 
lengths b = = b B , SCFT predicts a universal phase diagram [T8J [TS] that depends on 
only \N and the volume fraction Ja of the A block. 

We focus in this contribution on the behavior of the structure factor S(q) in simula- 
tions of symmetric diblock copolymer melts. S(q) is a correlation function for composition 
fluctuations, defined by 

S(q) = jdr (^(r)^(0))e^ r , (1) 

where ^(r) = [6ca{t) — 5ce(r)]/2 is a composition field. Here, <5cj(r) is a deviation of the 
number concentration of i monomers from the spatial average value = c/2 (c: total average 
monomer concentration). 

It is possible to use the SCF approximation to calculate S(q) , by employing the relation- 
ship between correlation functions and functional derivatives of F . This yields the so-called 
random phase approximation (RPA) [18J. The RPA yields a prediction for the quantity 
cNS~ 1 (q) in symmetric diblock copolymer melts of the form 

cNS-\q)=F(qR g0 )-2x e N (2) 

where F(x) is a known analytic function, R g0 = \J Nb 2 / 6 is the unperturbed radius of 
gyration, and Xe is an effective interaction parameter. In the most general form of SCFT, 
the effective interaction parameter that appears in the RPA is given by a derivative Xe = 
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(2kBT)~ 1 d 2 f ex (4>A) /d(j) 2 A of the SCF excess free energy per monomer. In the special case 
of symmetric diblock copolymer melts, the RPA thus predicts that S(q) depends on one 
non-dimensionalized wavenumber qR g o and a single thermodynamic parameter x e iV. 

In any theory that is based on the standard model, physical predictions must be invariant 
under changes in how we define " monomer" . The standard model is based upon a descrip- 
tion of polymer conformations as idealized fractals (Wiener processes), with no chemical 
microstructure. As such, the definition of a monomer, or the number N of monomers per 
chain, is arbitrary. Consequently, all SCF and RPA predictions can be expressed as relation- 
ships among dimensionless quantities that are invariant in this sense. Among such invariants 
are the product XeN, which is a measure of excess free energy per chain for transferring 
an A chain of N monomers into a B-rich environment, the invariant scattering intensity 
S(q)/(Nc), and quantities such as qR g o in which a length scale is non-dimensionalized by 
comparing it to a measure of the coil size. This invariance principle also, however, admits 
one parameter that does not appear in SCFT predictions (except as a trivial prefactor in 
expressions for total free energy), which is a measure of dimensionless concentration. In an 
incompressible one-component liquid of chains of length N, with a monomer concentration 
c and a polymer concentration c/N, we may define an invariant concentration 

C = cR 3 eQ /N = N 1 ' 2 ^ (3) 

which is the number of chains in a volume R^ , where R e o = vNb is a random- walk approx- 
imation for the root-mean-squared end-to-end length. Alternatively, following Fredrickson 
and Helfand [8], we may define an "invariant degree of polymerization" 

N = C 2 = N(cb 3 ) 2 . (4) 

We use R e o to define C and N, rather than the radius of gyration or some other measure of 
coil size, simply for consistency with the definition of N used in previous work [8j [161 02] • 
This choice is arbitrary. Notably, SCF predictions depend non-trivially upon all of the 
dimensionless parameters that are allowed by the invariance principle except the parameter 
C or N. 

The predicted absence of any dependence on N is a peculiarity of the SCF approximation. 
A series of closely related one-loop theories [H [121 H31 HS1 HZ] of the standard model, which 
go beyond the SCF approximation, all yield predictions that depend upon N as well as the 
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parameters that appear in SCF predictions. Moreover, these theories all predict results that 
reduce to corresponding SCF or RPA predictions in the limit N — > oo. These studies thus 
suggest that the SCF approximation is exact in the limit of infinitely strongly overlapping 
chains, but not finite chains. For symmetric diblock copolymer melts, the renormalized 
one-loop (ROL) theory predicts a structure factor that is given by a universal function 

cNS-\q) = H(qR g0 , X eN,N) (5) 

of qRgo, XeN, and N, where R g0 = ^N/Qb. The renormalized one- loop theory [TJl dSl 
[T7] yields an explicit prediction in which the function H is given as the sum of an RPA 
prediction (the RHS of Eq. ^ with renormalized values of Xe and b, plus a correction that 

1/2 — 1/2 

is proportional to iV . One-loop theories of fluctuations predict 0(N ) corrections to 
SCF predictions appear in many other quantities, including the critical region for polymer 
blends and corrections to the random walk statistics. It appears that this one-loop correction 
is simply the first term in a renormalized loop expansion of corrections to RPA and SCF 
predictions that (if pursued further) would yield a Taylor expansion of H(qR g o, XeN, N) in 
powers of N , in which the ROL theory for S(q) retains only the first two terms [16] . We 
thus expect the existence of a scaling relation of the form given in Eq. ^ to have a wider 
range of validity than the one-loop approximation. 

Extensive simulations of diblock copolymer melts have been carried out previously using a 
variety of coarse-grained lattice and continuum models. The use of coarse-grained models for 
this purpose presumes a high degree of universality, both among different simulation models 
and between coarse-grained simulations and experiments. Because none of these coarse- 
grained models provide an accurate representation of any real polymer liquid, it would be 
pointless to study them unless one believes that the equivalent behavior can be obtained 
from corresponding states of models or physical systems that differ in many details. The goal 
of the present contribution is to test the extent to which data for S(q) from several different 
simulation models of diblock copolymers is consistent with the existence of a universal 
function of the form given in Eq. 

The scaling hypothesis of Eq. ^ simply postulates that the invariant scattering intensity 
S(q)/cN exhibits a universal dependence upon all of the invariant dimensionless parameters 
that appear in the standard model. As such, a test of the adequacy of Eq. ^ is also a test 
of the adequacy of the standard model as a generic description of models and experimental 
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systems that differ in many microscopic details. 



III. SIMULATION METHODOLOGY 

We have studied three coarse-grained models. The first model is a lattice bond fluctuation 
model, which has been described in Refs. [22H21]- The second is a bead-spring model with a 
steeply repulsive Weeks- Chandler- Anderson (WCA) pair repulsion, similar to that used by 
Kremer and Grest [25]. The third is a bead-spring model with a much softer pair interaction 
that allows monomers to overlap [26T - T28] . We will refer to these models in what follows as 
models L ("lattice"), H ("hard"), and S ("soft"), respectively. 

Model L (the bond fluctuation model) is based on a cubic lattice, with a lattice spacing 
a, in which each monomer occupies (i.e., excludes other monomers from) the eight lattice 
sites of a cube of volume a 3 . The bond between consecutive monomers along a chain is 
allowed to have any of 108 vectors with lengths up to \/l0a, excluding the vector of length 
v8cr in the (2,2,0) direction to avoid bond crossing. Contacts between monomers of type i 
and j are assigned an energy cost e^, with €aa = £bb = —£ab = &/2. Simulations of this 
model used a monomer concentration ca 3 = 1/16, corresponding to half of the maximum 
allowed filling. Model L was simulated using lattice Monte Carlo, using a combination of 
local moves and of chain flip moves that interchange the A and B blocks of a chain. 

Model H used here is identical to that used in Ref. [IT] . In this model, all monomers 
interact via a purely repulsive Weeks-Chandler-Anderson (WCA) [29J pair potential, for 
which 

fy(r) = At^a/x) 12 - {a/xf + 1/4] (6) 

for i,j = A,B, with Vy(r) = beyond a cut-off r c = 2 1 / 6 cr, with €aa = ^bb = ksT and 
(-ab = £aa + a.. Bonded monomers interact via a harmonic potential fbond(^) = k{t — l) 2 /2, 
with k = 400/cbT/ct 2 and I = a. All simulations of this model were carried out in NVT 
ensemble in a L x L x L periodic cubic simulation cell, using a monomer concentration 
c = 0.7cr~ 3 . Data for this model was shown previously in Ref. [T7]. The model was 
simulated using a replica-exchange hybrid MC algorithm, in which short MD runs are used 
as proposed MC moves, among other MC moves, as described in Ref. [T7j . 
Model S uses a soft repulsive interaction 

V ij {r) = e ij [{r/a)-l) 2 (7) 
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cut off at r c = a. This was originally introduced by Groot and coworkers [26] in DPD 
simulations. We do not use a DPD thermostat here, instead we use a Nose-Hoover thermo- 
stat. Bonded particles in this model experience an additional harmonic bond potential with 
zero rest length and a spring constant k = 3.406A;bT/(T 2 . Simulations of the model used a 
monomer concentration c = 3.0a~ 3 , more than four times that of model H, thus creating 
substantial overlap between monomers. Model S was also simulated using a replica exchange 
hybrid MC algorithm (Simpatico [30]), using hybrid MC and chain flips moves, in which the 
HOOMD-blue code [3T] code was used as an engine to accelerate the hybrid MD moves on 
graphics processing units. 

In each of these three models, the pair interaction between monomers of types i and j 
is proportional to a parameter e^, with €aa = £bb- For each model, we define a parameter 
Oi = £ab — £aa, which we use as a control parameter to vary Xe- 

The structure factor in a diblock copolymer melt exhibits a peak at a wavenumber q* with 
a peak intensity S(q*) that increases with increasing degree of thermodynamic repulsion 
between A and B monomers, or increasing \ e . An example of this behavior, from our 
simulations of model S, is shown in Fig. 1. In experiment, Xe can be controlled by varying 
temperature. In the simulations presented here, we vary Xe by varying the control parameter 
a , while keeping the temperature and all other parameters of the Hamiltonian fixed. For 
each model and chain length, we have conducted simulations over a range of a > 0, up to 
values near those at which an ordered phase forms. We have determined values for the peak 
wavenumber q* and peak intensity S(q*) for each simulation by fitting S(q*) to a smooth 
function of q. An example of the quality of the fit is shown in Fig. 1. 

Our analysis requires a value of the statistical segment length b, which is used to define 
the parameters R g o = \JN/Qb and N = N(cb 3 ) 2 . A value of b was determined for each model 
by the procedure used in Refs. [151 E] an d in earlier work [32] by Wittmer and coworkers. 
For each model, we measure an A-dependent apparent segment length b 2 (N) = 6R 2 /N 
in a homopolymer reference state (with a = 0), for several different values of the chain 
length N. The chain length dependence of b(N) is very well described by the renormalized 
one-loop theory for mono-disperse homopolymer melts. We define the parameter b as the 
limit b 2 = liniAr_ s . 00 b 2 (N). This definition is required for consistency with the renormalized 
one-loop theory, which predicts random walk behavior only in this limit. This yields values 
of b = 3.244a for model L, b = 1.404a for model H and b = 1.088er for model S. 

7 




q^go 

FIG. 1. Reduced structure factor S(q)/cN as a function of the dimensionless wave number qR g o 
of a symmetric diblock copolymer melt for model S for N = 32 (N = 480), for a range of values 
a = £AB — £AA = 0, . . . , 1.35 (from bottom to top). 

The scaling hypothesis predicts that we should obtain equivalent behavior from different 
simulation models, over a range of values of a, from simulations of systems with equal values 
of N. To test this, we have designed pairs of simulations of different models in which the chain 
lengths have been adjusted to give matched values of iV = N(cb 3 ) 2 . Because different models 
have different values of the dimensionless parameter cb 3 , this requires that we compare 
simulations of different models with different values of N, the number of monomers per 
chain. To facilitate the design of simulations of models H and S with matched values of N, 
the value of the spring constant parameter k used for model S was chosen so as to yield a 
value of (cb 3 ) 2 = 14.9 that is exactly four times the value of this parameter obtained for 
model H. Simulations of model S with chains of iV monomers thus have the same value of 
iV as that obtained in simulations of model H using chains with 4N monomers. Below, we 
thus compare simulations of model S with chains of length N = 16 and N = 32 to model H 
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simulations of much longer chains of length N = 64 and N = 128. 



IV. TESTING UNIVERSALITY 

To test the scaling hypothesis of Eq. (J5]), we compare results for the non-dimensional peak 
wavenumber q*R g o and invariant peak intensity cNS" 1 ^*) for the three simulation models 
described above. Eq. (|5| implies that these quantities should both be universal functions 

q*R g0 = Q*(XeN,N) (8) 
cNS~\q*) = H*( X eN,N) (9) 

of the two variables XeN and N, where R g o = y N/6b. The RPA predicts a constant value 
of q*R g o = 1.95, independent of XeN, and a value for H* that varies linearly with XeN, 
independent of N. As part of our hypothesis, we interpret the microscopic parameters Xe and 
b for each model as "renormalized" parameters that can exhibit an arbitrary dependence on 
temperature, monomer, concentration and on the details of the monomer-scale interactions, 
but that are assumed to be independent of chain length. Specifically, we assume that the 
value of Xe for each of our simulation models is some unknown, generally nonlinear function 
Xe{oi) of the control parameter a, and that the function Xe( a ) is different for different models. 

The scaling hypothesis predicts that we should obtain equivalent behavior in correspond- 
ing states of different models. Corresponding states of symmetric diblock copolymer melts 
are defined as states with equal values of the N and of the product XeN. A useful alter- 
native definition corresponding states may be obtained by defining an invariant interaction 
parameter Xei D Y requiring that x e N = XeN. This definition yields 

X e = X e/(cb 3 ) 2 . (10) 

Corresponding states can then also be described as states with equal value of x e an d N, 
which is equivalent to requiring equal values of x e iV and N. We assume that % e , like Xe, 
is a model-dependent function of a, but is independent of N. This definition is needed in 
subsection |IV C , where we define corresponding values of a in different models by requiring 



that they yield equal values for x e - 

In what follows, we always compare results of pairs of simulations of different models in 
which the values of N have been chosen to give equal of N. After creating simulations with 



matched values of N, the remaining conceptual difficulty that must be overcome to test 
Eqs. ^ and ^ is the the fact that the equations postulate a dependence on a parameter 
Xe{oi) that is an unknown function of the control parameter a, and that cannot be directly 
measured. Two strategies are used in what follows to overcome this difficulty. The first, 
which is described in subsection A, is to test the universality of the relationship between 
two measurable quantities, the peak wavenumber and peak intensity. The second, which is 
described in subection C, involves an attempt to collapse the data from different models, 
over a range of chain lengths, by constructing an explicit mapping between corresponding 
values of a in different models. 

A. Peak wavenumber vs. peak intensity 

As our first test of universality, we construct parametric plots of the relationship between 
observables q*R g o and cAS ,_1 (g*) for pairs of models with matched values of N, over a range 
of values of the control parameter a for each model. Eq. ^ implies that systems with equal 
values of both N and cAS' _1 (g*) must be in corresponding states, with equal values of XeN. 
Systems with equal values of N and cNS~ l (q*) should thus also exhibit equal values for the 
peak wavenumber q*R g o- This implies that a parametric plot of q*R g o vs. cAS ,_1 (g*) should 
yield a curve whose form depends on N, but that data from different simulation models with 
identical values of N should collapse onto the same master curve. 

Figures [2] and [3] show plots of q*R g o vs. a parameter x* a N that is directly related to the 
scattering intensity, for several pairs of simulations of different models with matched values 
of N . The quantity x* is an "apparent" interaction parameter, also used in Ref. [17], that 
we define to be the value of x that would be inferred by fitting the measured peak intensity 
to the RPA prediction, by requiring that 

cNS-\q*)=2[( X N) c -xlN] (11) 

where {xN) c = F(q*R gQ )/2 = 10.495 [IS] is the critical value of XeA for symmetric diblock 
copolymers, at which the RPA predicts a divergence in S(q*). Plotting qR g o vs. x* a N is 
thus equivalent to plotting of qR g o vs. cNS~ 1 (q*), but show the data in a form in which 
peak intensity increases from left to right, and in which x* a N = 10.495 corresponds to a 
divergence of S(q*). The RPA also predicts a constant value of q*R g o = 1.95 that depends 
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on neither x* a N nor N, which is shown by a horizontal dashed line in these figures. 

Figure [2] shows simulation results for q*R g o vs. x* a N for models L and H, for pairs of 
simulations in which the values of N have been chosen to give nearly matched values of N 
(to the closest integer value of N). The upper panel shows a comparison of simulations with 
N ~ 75 and the lower panel shows a pair of simulations with N ~ 145. When we compare 
data from simulations using matched values of N, the level of agreement between the lattice 
and continuum model is remarkable. 

Figure [3] shows data for q* /q vs. Xa,N for the H (circles) and S (triangles) models, for two 
pairs of simulations with exactly matched values N = 240 (upper panel) and iV=480 (lower 
panel). In each of these pairs of simulations, matched values of N are obtained by using 
chains in model H that have 4 times as many monomers as those used in the corresponding 
model S, because of the much higher monomer concentration in the latter model. Here, 
data for each model is shown for two different values of the simulation box length L, to 
verify the absence of finite size effects. Agreement between the different models is again 
nearly perfect. This collapse of the data for models with such widely disparate chain lengths 
provides a particularly clear confirmation that the form of this function depends only on 
overlap parameter N, and not on the actual number of monomers. 

B. Failure of Perturbation Theory 

In previous work [T7j, we compared results for the invariant peak intensity for model H 
to several different theories by using perturbation theory to estimate Xe(°d- I n Ref- [IS], it 
was shown how Xe{ a ) could be calculated to first-order in a Taylor expansion in a, of the 
form 

Xe (a) = z°°a/k B T + 0(a 2 ) (12) 

in which the value of the coefficient z°°, known as the effective coordination number, can 
be obtained from simulations of a reference homopolymer (a = 0) melt. This analysis 
generalized and provided a more rigorous statement of an earlier proposal by one of us 
(M.M.) and Kurt Binder [231 E3] that Xe be estimated for the lattice model by taking 
Xe{oi) = z(N)a/kBT, where z(N) was the number of inter-molecular neighbors of each 
site in a lattice simulation of a reference homopolymer melt containing chains of length 
N. The corresponding definition of z(N) for a continuum model is also a measure of the 
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FIG. 2. Nondimensional peak scattering location q*R g o, vs. apparent interaction parameter x* a N 
for simulations of model L (the lattice bond fluctuation model) and model H (a bead spring model 
with harsh repulsion). Data is shown for two pairs of simulations with nearly matched invariant 
degrees of polymerization N ~ 74 (upper panel) and N ~ 145 (lower panel). The RPA prediction 
q*R g o = 1-95 is shown as a horizontal dashed line. 
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FIG. 3. Normalized peak scattering location q*R g o vs. x* a N for model H (open symbols) and model 
S (filled symbols). Data is shown for two pairs of simulations with matched values of N = 240 
(upper panel) and N = 480 (lower panel). For each model, data is shown for two different values 
of the simulation box dimension L, using circles for the smaller box and squares for the larger, as 
indicated in the legend. The RPA prediction is shown by a dashed horizontal line. Data for model 
H are shown over ranges a = - 1.2 for N = 64 and a = - 0.4 for N = 128. Data for model 
model S are shown over ranges a = - 6.4 for W = 16 and a = - 4.7 for N = 32. 
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FIG. 4. Normalized inverse peak intensity (l/2)cNS~ 1 (q*) for the lattice bond fluctuation (L) 
model (lines and open squares) and the continuum 'hard' (H) model (solid circles), for three pairs 
of simulations with nearly matched values of N = 75, 145 and 292. Data for both models is plotted 
here vs. xi^^ = z°°aN /ksT ', where = 0.2965 for model H and z£° = 4.2 for model L. 



number of inter-molecular neighbors of each monomer, which is expressed in this case as an 
integral involving the inter-molecular radial distribution function. Ref. [15] clarified that, 
in detailed comparisons with theory, the SCF parameter \ e must be defined using a value 
z°° that is denned by extrapolating z(N) to iV = oo. Following the prescription outlined in 
that reference, we have obtained values of z°° = 4.2 for the lattice model, z°° = 0.2965 for 
model H, and z°° = 0.237 for model S. 

Comparisons of data for cNS^ 1 (q*) vs. XeN to the (ROL) theory, using this estimate 
for Xe, showed good agreement at low values of a, but exhibited systematic deviations at 
higher values of a [17]. Because large deviations occurred only at large values of a, they 
were tentatively ascribed to a breakdown in the first-order expansion for Xe{ot)- In that 
study, however, we could not cleanly separate failures of the Taylor expansion of Xe(«) from 
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failures of the ROL theory. 

A more basic test of the adequacy of a first-order perturbation theory for Xe( a ) may be 
obtained by using this estimate for Xe(oc) to compare two different models, using pairs of 
simulations with matched values of N. Fig. [4] shows a comparison of data for ciVS 1-1 ^*) 
vs. Xe( a ) f° r three pairs of simulations of models L and H, plotting using the estimate 
z°°aN/kBT, for values of N ~ 75, 145 and 245. When plotted in this way, the results from 
the two models disagree rather badly, and agree only at low values of XeN. This indicates a 
severe failure of either the scaling hypothesis or of the first-order approximation for Xe( a )- 
In the next subsection, we show that it is the latter. 

C. Mapping between corresponding states 

To test Eq. ^ without making any assumptions about how Xe depends on a, we have 
tested whether it is possible to collapse the data for different models by constructing an 
explicit nonlinear mapping between corresponding values of the control parameter a in 
different models. 

Consider a comparison of two models, 1 and 2, with control parameters ai and «2- We 
assume that x e * s given by different functions of a in the two models, but that there exists 
a mapping between corresponding values of a of the two models. We define corresponding 
values of a so as to give equal values for x e - Hence, corresponding values of a and N 
yield corresponding states. We assume the existence of a mapping Q{a,2) such that using 
a.\ = Q(ac2) in model 1 yields the same value for x e an d as that obtained by using «2 
in model 2. Given data for the invariant inverse peak intensity ciVS 1-1 ^*) from a pair of 
simulations of these two models with matched values of N, over a range of values of a, it 
should be possible to collapse this data by plotting the results from model 1 vs. a,\ and 
plotting the results from model 2 vs. the corresponding quantity Q(a 2 ). 

Because we do not know the mapping function Q a priori, we must test universality 
by searching for a a mapping function that yields the best possible collapse of results for 
cNS~ l (q*) from pairs of simulations with matched values of N. The required fitting proce- 
dure would be meaningless if applied to only a single pair of simulations, with a single value 
of N: If cNS~ 1 (q*) is a monotonic function of a, there will always exist a mapping that 
collapses two such curves, since it will always be possible to identify pairs of values of a that 
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yield the same invariant intensity. The procedure provides a meaningful test, however, if it 
is applied to several such matched pairs of simulations of the same two models, with a range 
of values N, by attempting to simultaneously collapse data for several pairs of simulations 
with different values of N using a single mapping Q between corresponding values of a. 

The graphical construction described above requires that we plot the invariant intensity 
vs. the control parameter a of one of the two models. This choice of coordinate is clearly 
not unique. If we can find a mapping Q that collapses the data in this representation, then 
we can also plot the data for model 1 vs. any function f(a\) and then plot data for model 
2 as a function of a quantity f(Q(o(2)), without changing the quality of the collapse. 

In practice, we have chosen to plot the data for ciVS 1-1 ^*) for one model in each pair 
(model 1) as a function of the quantity NixiP = zfaNi/kgT used in Fig. [4] to test the first 
order perturbation theory for x&\ an d to plot the data from the other model (model 2) as a 
function of a quantity A^ 2 x™ ap (a2); i n which X™ ap (a2) is a polynomial whose coefficients are 
fitted such as to collapse the data. The quantity N 2 xT ap yields the value of NiXe^ in model 
1 that corresponds to the same value of x e) or XeN, as that obtained by using a 2 in model 
2, where N\ and N 2 denote values of N in simulations with matched values of N. To first 
order in a, x™ ap (o0 is given by x™ ap ( a ) — z^a/ksl \ where z°° is the effective coordination 
number for model 2. In both of the comparisons shown below, we have fit xT ap ( a z) using a 
third polynomial of the form x™ ap ( a ) = z°°a + Bo? + Ca 3 in which the coefficients B and 
C are adjusted to collapse the data. 

Figure [5] shows the results of such an attempt to collapse the peak intensity data for 
models L (open squares) and H (open circles). The values of N and N are the same as those 
used in Fig. |4j Here, data for the continuum model is plotted as a function of Xe(cn)N while 
using the first-order approximation Xe(&) = z^a/ksT for Xe{oi). Data for the lattice model 
is plotted vs. iVx™ ap (a), using a third order polynomial fit for x™ ap (a), as described above. 
In contrast to the comparison shown in Fig. [4| the collapse is nearly perfect. This result 
confirms that the results of these two models can be adequately described by a single scaling 
function, but that it is essential to allow for a nonlinear relationship between corresponding 
values of a in different models. 

An analogous comparison of results for models H and S is shown in Fig. |6j for two pairs 
of simulations with matched values of iV = 240 and iV = 480, as in Fig. [3] Here, the inset 
shows a comparison in which data from models is plotted using the first order perturbation 
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FIG. 5. Collapse of the invariant inverse peak intensity (l/2)cNS~ 1 (q*) for model L (the lattice 
model, lines and open squares) and model H (continuum model, solid circles), for three pairs of 
simulations with nearly matched values of N = 75, 145 and 292. The data are the same as in 
Fig. [4J Data for model H are plotted here vs. xi 1 ^ = z°°aN /ksT, where z°° = 0.2965. Data for 
model L are plotted vs. Xe ap (ct)N , using a cubic polynomial Xe &P = z^a + 4.94a 2 + 21.46a 3 in 
which the coefficients of the quadratic and cubic terms have been adjusted to collapse the data. 



theory estimates for Xe(oc), and the main plot shows an attempt to collapse the data by 
allowing for a nonlinear mapping. In the main plot, data for model H is plotted using the 
first order perturbation theory for Xe( a ), an d data for model S is plotted vs. iVx™ ap (cv), 
were xT ap ( a ) ^ s a third order polynomial with two coefficients that are chosen to collapse 
the data. The difference between the models in the inset is less dramatic than for models H 
and L. Allowing for a nonlinear mapping again allows us to obtain nearly perfect collapse 
of data from two models. 
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FIG. 6. Collapse of the invariant inverse structure factor cNS^ 1 (q*) for model H ('hard') and 
model S ('soft'), for the same set of simulations as those shown in Fig. [3j The inset shows results 
for both models plotted using the first-order approximation x^e\ a ) = z 00 ajksT . In the main plot, 
data from model H is plotted vs. the first-order approximation for Xe\<x)N and data for model S is 
plotted vs. x^ lap (a)iV, where x™ ap (o;) is a third order polynomial x™ ap (a) = Zg°o; + 0.01a 2 + 0.03a 3 
in which two coefficients have been adjusted so as to collapse the data. 



V. CONCLUSION 



Theories that are based on the highly idealized standard model of polymer liquids implic- 
itly assume the validity of a principle of corresponding states: They presume that equivalent 
behavior should be obtained from simulations and experimental systems that can differ in 
many microscopic details, as long as we compare systems with equal values for all of the 
dimensionless parameters that appear in the standard model. The work presented here is 
motivated by the observation that this is a testable scaling hypothesis, and which can be 
tested without reference to any quantitative theory, by comparing results of different simu- 
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lation models. In this paper, we have tested this scaling hypothesis by comparing results for 
the structure function S(q) from three substantially different simulation models of symmetric 
diblock polymer melts. The results strongly support the scaling hypothesis. 

The main conceptual difficulty that had to be overcome was the fact that the hypothesis 
allows all of the quantities of interest to exhibit an unknown dependence on a phenomenolog- 
ical interaction parameter that is an unknown, model-dependent function of the parameters 
that one can control in a simulation. We devised two tests of universality that did not 
require knowledge of this dependence. The first of these involved a comparison of the rela- 
tionship between two observable quantities, the peak wavenumber and the peak intensity. 
The second involved a test of whether the data for peak intensity from different models could 
be collapsed by searching for an iV-independent mapping between corresponding values of 
the control parameters in different models. We believe that analogous methods of analyzing 
data may be useful for analysis of experimental data. In order to compare to experiments, 
however, it may be important useful to allow for the effects of various types of asymmetry 
that are not present in these simulations. 

Our results strongly support the scaling hypothesis. Results for the peak position and 
peak intensity from three substantially different models were found to exhibit almost perfect 
universality. For comparisons involving model S, this was found even for surprisingly short 
chains, with as few as 16 monomers. Notably, the level of agreement between different models 
for the relationship between the peak wavenumber and peak intensity was substantially 
better than the already rather good agreement obtained between results of model H and the 
ROL theory [17] . This is consistent with the idea that the accuracy of the scaling hypothesis 
is a measure of the accuracy of the standard model itself, which should have a wider range 
of validity than the one-loop approximation for this model. 

We emphasize that the goal of this paper was not to test the accuracy of a particular 
theory for S(q) or other equilibrium properties. Nothing in the analysis presented here 
referred to a specific theory. Instead, we sought to test whether it was possible for any coarse- 
grained theory based on the standard model to describe the data from different simulation 
models. One consequence of this approach is that our analysis did not yield estimates 
the interaction parameter \ e for any of the models that we have studied. The only way 
to obtain an estimate of Xe for a particular model or experimental system is to compare 
simulation or experimental results to a quantitative theory. The quality of the resulting 
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estimate can, of course, be no better than the quality of the theory upon which it is based. 
The analysis presented here shows only that it is possible to construct a theory (or simply 
a set of empirical correlations) that can describe a variety of simulation models. 

The scaling hypothesis that we discuss here can be easily generalized to ordered phases 
of block copolymer melts, and to other types of dense polymer liquids, such as polymer 
mixtures. In future work, we plan to extend this analysis so as to test the universality of 
behavior in the vicinity of the order-disorder transition of diblock copolymer melts, and in 
the ordered phases. 

The hypothesis that we tested here is intentionally very simple, in that it allows physical 
properties to depend only on the variables that appear explicitly in the standard model. It 
neglects a variety of phenomena that could lead to non-universal deviations from scaling in 
more realistic models and in experimental systems. Among these are the possible existence 
of a square-gradient corrections to the local excess free energy, leading to a non-universal 
q 2 contribution to Xe, and the possibility of special corrections to the excess free energy 
arising from chain ends and junctions, as discussed in Ref. [12] . The fact that the scaling 
hypothesis holds to very high accuracy in the models studied here implies, however, that 
such non-universal corrections must be very small in this particular set of models. This is 
encouraging, but it remains to be seen how large may be the deviations from scaling in other 
simulation models, and in experiments. The type of analysis presented here is, however, the 
natural starting point for quantifying any such deviations. 
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